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Bounds on the Higgs mass from the Tevatron and LHC are determined using exclusive jet bins 
to maximize sensitivity. Scale variation in exclusive fixed-order predictions underestimates the 
perturbative uncertainty for these cross sections, due to cancellations between the perturbative 
corrections leading to large K factors and those that induce logarithmic sensitivity to the jet-bin 
boundary. To account for this, we propose that scale variation in the fixed-order calculations should 
be used to determine theory uncertainties for inclusive jet cross sections, whose differences yield 
exclusive jet cross sections. This yields a theory correlation matrix for the jet bins such that the 
additional uncertainty from large logarithms due to the jet boundary cancels when neighboring 
bins are added. This procedure is tested for _ff -I- 0, 1 jets, WW + jets, and W + 0,1,2 jets, 
and found to be generally applicable. For a case where the higher-order resummation of the jet 
boundary corrections is known, we show that this procedure yields fixed-order uncertainties which 
are theoretically consistent with those obtained in the resummed calculation. 



I. INTRODUCTION 

In the search for the Higgs boson at the Tevatron and 
the Large Hadron Collider (LHC), the data are divided 
into exclusive jet bins. This is done because the back- 
ground composition depends on the number of jets in the 
final state, and the overall sensitivity can be increased 
significantly by optimizing the analysis for Higgs -I- 0, 1, 
and 2 jet signals. The primary example is the H — >■ WW* 
decay channel, which dominates the current Tevatron ex- 
clusion limits around ttih — 2mw [ll 12] , and is one of the 
important channels for tuh ^ 130 GeV being pursued at 
the LHC 0, |i|. The importance of the Higgs + 1 jet 
channel in H ^i- tt and H -^ WW* was demonstrated 
explicitly in Refs. [5|, |6|. Similarly, for H -^ 77, which 
plays an important role for ttih ^ 130 GeV, the search 
sensitivity can be improved by optimizing the analysis 
for different jet bins [7]. 

Since the measurements are performed in each jet bin, 
the perturbative uncertainties in the theoretical predic- 
tions must also be evaluated separately for each jet mul- 
tiplicity Q. Furthermore, to combine the results in the 
end, the correlations between the theoretical uncertain- 
ties in the different jet bins as well as in the total cross 
section have to be understood and taken into account. 

In the winter 2011 Tevatron analyses of gg — > iJ — > 
WW* [2], the perturbative uncertainties in the signal 
cross section are evaluated using common scale variation 
for the exclusive jet bins, which yields jSj] 

^^^ = 66.5% X (+^^:) + 28.6% X (+^4%) 

+ 4.9%x(+I«^:) = (;-J:). (1) 

The three terms are the contributions from the 0, 1, and 
(> 2)-jet bins with their relative scale uncertainties in 
brackets. By using a common scale variation the uncer- 
tainties are effectively 100% correlated and are added lin- 



early, such that the ±14% scale uncertainty in the total 
cross section is reproduced. 

For the 0-jet bin, which is the most sensitive search 
channel in H —^ WW*, one applies a strong veto on 
additional jets. It is often argued that with the jet veto 
the perturbative uncertainties improve, yielding scale un- 
certainties from fixed-order perturbation theory that are 
smaller than those in the total cross section, as seen in 
Eq. ([1]). This apparent improvement arises from can- 
cellations between two sources, large corrections to the 
total cross section (large K factors) and the large correc- 
tions from logarithmic dependence on the jet veto. Since 
the improvement arises from a cancellation between two 
large and predominantly independent perturbative series, 
it must be assessed carefully. 

We propose a simple procedure to estimate more real- 
istic perturbative uncertainties for exclusive jet bins from 
fixed-order perturbation theory. The method is designed 
for processes with large K factors or large perturbative 
corrections in inclusive cross sections and takes into ac- 
count the structure of the various perturbative series. As 
we will see, it can also be applied in general. The essential 
idea is to first independently determine the uncertainties 
in the inclusive A^-jet cross sections cr>Ar, and then use 
them to compute the uncertainty in the exclusive -/V-jet 
cross section a^ from the difference 



CTat = cr>jv — cr>jv+l 



(2) 



To a first approximation the perturbative series for a>N 
can be considered unrelated for different N. For instance, 
their series start at different orders in as, and there is 
a priori no direct relation between the modifications to 
the series caused by the jet cuts that define these two 
inclusive samples. Therefore, as explained in detail in 
Sec. [ni we can work in the limit where the fixed-order 
perturbative uncertainties in the cr>Ar's can be taken as 
uncorrelated, leading to 



A2 



A2 



A2 



(3) 



The uncertainty in the exclusive cross section is larger 
than that in the corresponding inclusive one, which ac- 
counts for its more complicated perturbative structure. 
Equation ^ also leads to an anticorrelation between the 
cross sections in neighboring jet bins. When neighboring 
bins are added the sensitivity to the boundary between 
them cancels and the uncertainty reduces accordingly. 

For example, for the 0-jet bin m H —>■ WW* discussed 
above, we have ctq = crtotai — <y>i- Here, tT>i contains 
double logarithms of the jet pr cut, whereas CTtotai does 
not involve any jet definition, so their perturbative series 
can be considered largely independent. Therefore, tak- 
ing their perturbative uncertainties Atotai and A>i as 
uncorrelated, the covariance matrix for {crojO'>i} is^ 



A^ 

"total 



A^>i 



-A2 



All. 



Using this matrix to compute the uncertainty in cro + cr>i 
reproduces Atotai as it should. 

We should mention that we are only discussing here 
the uncertainties due to unknown higher-order pertur- 
bative corrections, which are commonly estimated using 
scale variations. We do not discuss parametric uncer- 
tainties, such as parton distribution function (PDF) and 
as uncertainties, which have been extensively discussed, 
recently for example in Refs. [9l-[l6J. 

In the next section we present the arguments leading 
to our proposal for evaluating the perturbative uncertain- 
ties for exclusive jet bins, and discuss the structure of the 
perturbative series. In Sec. IIIIl we apply our method to 
a variety of processes. We start in Sees. IIII Al and IIII Bl 
with discussion and numerical results for gg -^ H + 
jets and gg —> H + I jets. In Sec. IIII CI we consider 
pp -^ WW+0 jets, which is an important background for 
Higgs production. In Sees. IIII Dl IIII El and IIII Fl we con- 
sider W + 0,1,2 jets, which are important backgrounds 
for missing-energy searches. In Sec. lIVI we consider again 
gg ^ H + jets and test our method for the fixed-order 
uncertainties against a case where the resummation of 
the large logarithms induced by the jet binning is known 
to next-to-next-to-leading logarithmic (NNLL) accuracy. 
We conclude in Sec. [V] In the Appendix, we give expres- 
sions for the uncertainties and correlations for the case 
where one considers 0, 1, and (> 2)-jet bins as in Eq. ([T]). 



II. JET BIN UNCERTAINTIES 

To examine in more detail the modification of the per- 
turbative series that takes place for exclusive jet bins, 



we will consider the example of the 0-jet bin and (> 1)- 
jet bin. The total cross section, crtotai, is divided into a 
0-jet exclusive cross section, ao{p^^^), and the (> l)-jet 
inclusive cross section, cr>i(p'^"*). 



Ctotal 



dp-T- 
ap 



^y Jp' 

fT0(p™*) + fT>l(p-*). 



dp—- 
dp 



(5) 



Here, p denotes the kinematic variable which is used to 
divide the cross section into jet bins. For most of our 
analysis we take p = p!^ , which for Eq. ([S]) is the largest 
Pt of any jet in the event. In this case, cro(PT'') only 
contains events with jets having pT < p't^, and o'>i(p™*) 
contains events with at least one jet with pT > p™' . 
In Eq. ([5]) both ctq and cr>i depend on the phase space 
(4) cut, p"^"', and by construction this dependence cancels 
in their sum. This means that the additional perturba- 
tive uncertainty induced by this cut, call it Acut, must 
be 100% anticorrelated between cro(p'^"*) and cr>i(p'^"'). 
That is, the contribution of Acut to the covariance matrix 
for {cto, (T>i} must be of the form 



Oc 



^A2 



(6) 



The questions then arc: (1) how can we estimate Acut, 
and (2) how is the overall perturbative uncertainty Atotai 
of (Ttotai related to the uncertainty for ao and a>i. 

To answer these questions, we discuss the perturba- 
tive structure of the cross sections in more detail. By 
restricting the cross section to the 0-jet region, one re- 
stricts the coUinear initial-state radiation from the collid- 
ing hard partons as well as the overall soft radiation in 
the event. This restriction on additional emissions leads 
to the appearance of Sudakov double logarithms of the 
form L"^ = In^ {p™*/Q) at each order in a perturbative 
expansion in the strong coupling constant as , where Q is 
the hard scale of the process. For Higgs production from 
gluon fusion, Q = ruH, and the leading double logarithms 
appearing at 0{as) are 



3a,, 
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ao(pr)=^B 1-— 21n^i:^ 
IT niH 



(7) 



where as is the Born (tree-level) cross section. 

The total cross section just depends on the hard scale 
Q, which means by choosing the scale ^ c:^ Q, the fixed- 
order expansion does not contain large logarithms and 
has the structure^ 



Ctotal 



(Jb[1 



0{al 



(8) 



^ Since these are theory uncertainties, there is no strict reason to 
combine them in a particular way. We add them in quadrature 
since this is the most convenient for discussing correlations and 
error propagation. 



^ These expressions for the perturbative series are schematic. They 
do not show the convolution with the parton distribution func- 
tions contained in ag, nor do they display ^ dependent loga- 
rithms. In particular, the single logarithms related to the PDF 
evolution are not displayed, since they are not the logarithms we 
are most interested in discussing. 



The coefficients of tfiis series can be large, corresponding 
to the well-known large K factors. For instance, the cross 
section for gg -^ H doubles from leading order to next- 
to-leading order (NLO) even though Ug ^ 0.1. As usual, 
varying the scale in as (and the PDFs) one obtains an 
estimate of the size of the missing higher-order terms in 
this series, corresponding to Atotai- 

The inclusive 1-jet cross section has the perturbative 
structure 

a>i(p™*)~aB[a.(L2 + L + l) (9) 

+ aliL^ + L^ + L2 + L + 1) + OialL"^)] , 

where the logarithms L = ln(p'^"'/(5) arise from cutting 
off the IR divergences in the real emission diagrams. For 
p'^"' <C Q the logarithms can get large enough to over- 
come the Us suppression. In the limit asL^ — 1, the 
fixed-order perturbative expansion breaks down and the 
logarithmic terms must be resummed to all orders in as 
to obtain a meaningful result. For typical experimental 
values of p*^"' fixed-order perturbation theory can still be 
considered, but the logarithms cause large corrections at 
each order and dominate the series. This means varying 
the scale in a^ in Eq. ^ directly tracks the size of the 
large logarithms and therefore allows one to get some es- 
timate of the size of missing higher-order terms caused 
by p'^"*, that correspond to Acut- Therefore, we can ap- 
proximate Acut = A>i, where A>i is obtained from the 
scale variation for cr>i. 

The exclusive 0-jet cross section is equal to the differ- 
ence between Eqs. ([5]) and ^, and so has the schematic 
structure 



is diagonal, 



ao(p^"*) ^aB{[l + as + a^s + 0(a^)] 
- [as{L^+L + l)+al{L'^- 
+ 0{alL')]}. 



L*+L^ 



L + l) 
(10) 



In this difference, the large positive corrections in crtotai 
partly cancel against the large negative logarithmic cor- 
rections. For example, at 0{as) there is a value of L 
for which the Ug terms in the schematic Eq. ([T0| cancel 
exactly, indicating that at this p'^"* the NLO cross sec- 
tion has vanishing scale dependence and is equal to the 
LO cross section, ao(p'^"*) = <tb- We will see this ef- 
fect explicitly in our examples below, using the complete 
perturbative expressions. We will find that this occurs 
for values of p'^"* in the experimentally relevant region. 
Because of this cancellation, a standard use of scale vari- 
ation in Eq. ([T0| does not actually probe the size of the 
logarithms, and thus is not suitable to estimate Acut- 

Since Acut and Atotai are by definition uncorrelated, 
by associating Acut = A>i we are effectively treating 
the perturbative series for crtotai and cr>i as independent 
with separate (uncorrelated) perturbative uncertainties. 
That is, considering {(Ttotai, o'>i}, the covariance matrix 



A2 

"total 





A2 



(11) 



This is consistent, since for small p'^"' the two series have 
very different structures. In particular, there is no reason 
to believe that the same cancellations in (Tq will persist 
at every order in perturbation theory at a given p"^"'. 

From Eq. ([TT|) it follows that the perturbative uncer- 
tainty in CTo(p'^"*) is given by A^^^^j + A>j, i.e., by sum- 
ming the inclusive cross section uncertainties in quadra- 
ture. It also follows that the complete covariance matrix 
for the three^ quantities {crtotai, ctq, cr>i} is 



^ a2 
"total 



A2 
"total 







c 



A 2 A 2 , A 2 _ A 2 

^total ^>1 ^^ ^total ^>1 



(12) 







-A2 



A^ 



where Atotai and A>i are considered uncorrelated and 
are evaluated by separately varying the scales in the 
fixed-order predictions for crtotai and cr>i(p'^"'), respec- 
tively. The A>i contributions in the lower right 2x2 
matrix for (j^y and cr>i are equivalent to Eq. ^ with 
Acut = A>i. Note that in this 2x2 space all of Atotai 
occurs in the uncertainty for uq. This is reasonable from 
the point of view that erg starts at the same order in as as 
ftotai and contains the same leading virtual corrections. 

The limit Acut = A>i which Eq. P^ is based on is 
of course not exact but an approximation. However, the 
preceding arguments show that it is a more reasonable 
starting point than using a common scale variation for 
the different jet bins. The latter usually results in the 
cross sections being 100% correlated, as in Eq. ([T]), and 
in particular does not account for the additional p'^"' in- 
duced uncertainties. In our numerical examples below, 
we will see that our method produces more sensible un- 
certainty estimates for fixed-order predictions. In SeclTVl 
we will compare the estimates from our method with 
those obtained by an explicit resummation in the jet- 
veto variable. This provides further evidence that our 
method gives consistent uncertainty estimates. Resum- 
mation provides a way for improving predictions for the 
central value of the cross section, together with better 
estimates of Acut and the structure of the theory corre- 
lation matrix, as discussed in Sec. IIVI 

It is straightforward to generalize the above discussion 
to jet bins with more jets. For the A'^-jet bin we replace 
CTtotai -^ cr>Ar, ctq -)■ cr^v, and cr>i -> cr>Ar+i, and take 
the appropriate ub- If the perturbative series for a>is! ex- 
hibits large as corrections, then the additional large loga- 
rithms present in cr>jv+i will again lead to cancellations 



^ The fact that only two of three are independent is reflected in 
the matrix, i.e. any 2x2 submatrix can be used to derive the 
full 3x3 matrix using the relation crtotai = co + cr>t. 



when we consider the difference ctat = (J>n — <y>N+i- 
Hence, A>7v+i will again give a better estimate for the 
Acut that arises from separating cr>jv into jet bins a^ and 
a>N+i- Another advantage of our procedure is that it is 
easily generalized to more than two jet bins by iteration. 
The case of three jet bins is given in the Appendix. 



III. EXAMPLE PROCESSES 

To elucidate the effect of pl^ vetoes on the fixed-order 
cross sections and demonstrate our method, we will now 
go through several explicit examples, considering in turn 
H + Q jets, H + 1 jet, WW + jets, and VF + 0, 1, and 
2 jets. All of our NLO pr spectra are obtained using 
the MCFM code [T7l-l20t. As our jet algorithm we use 
anti-fcr for the LHC results and a cone algorithm for the 
Tevatron results with R = 0.5 for both. 



A. Higgs + Jets 

In Higgs production via gluon fusion the cross section is 
known to next-to-next-to- leading order (NNLO) l2ll - l28l |. 
and exhibits large perturbative corrections. Consider the 
numerical results for the Higgs production cross section 
for ruH = 165 GeV, /x/ = /x^ = toh/2, and MSTW2008 
NNLO PDFs [p, for which as = as(mH/2) = 0.1189. 
Here one finds |30l - [33| 

atotai = (3.32 pb) [l + 9.5 a. + 35 a^ + 0[al)\ , (13) 

for the LHC at ^cm = 7TeV. Note that there is an 
a^ in the Born cross section, ub — 3.32 pb, but only 
the relative size of the corrections is important for our 
discussion. For the Tevatron the series is 

CTtotai = (0.15 pb) [1 + 9.0 as + 34 a^ + 0{al)\ . (14) 

In both cases the large K factors are clearly visible.** For 
the inclusive 1-jet cross section at the LHC one finds 

(Jyiip^T > 30 GeV, IV''*] < 3.0) 

= (3.32pb)[4.7a^ + 26a^ + 0(a^)] , 



<J>i{t^t > 25 GeV) 



(3.32 pb) [6.0 as + 32 a^ + 0{al)\ 



(15) 



Using instead /if = iir = mu the coefficients of the as and 
Og terms increase to 11 and 65 for the LHC and 12 and 74 for 
the Tevatron, respectively. The a^ coefficients for the Tevatron 
for example arise as 9.0 = 4.9 + 2.0 + 2.1 {/i = m,jj/2) and 
12.0 = 4.9+5.7+1.4 {/i = mn) where the three contributions arc 
respectively from the terms in the partonic cross section propor- 
tional to 5(1 — z), terms involving the plus functions [1/(1 — z)] + 
and [ln(l — 2;)/(l — z)]-!-, and the remaining terms that are nonsin- 
gular for z —^ 1. When separating these different terms we keep 
the overall 1/z factor in the convolution integral with measure 
dz/z. 



The first values correspond to the ATLAS and CMS ref- 
erence cuts, and the second to current ATLAS and CMS 
H —5- WW* analyses [1, |3| • Similarly, for the typical cuts 
used in H -^ WW* at the Tevatron 0, one finds 



cr>i{p'T >20GeV, iTyJ''*! < 2.5) 

= (0.15pb)[4.1a, -f 27a^ + 0(a^)] 



(16) 



In both Eqs. p3|) and (fTB|) one clearly sees the impact of 
the large logarithms on the perturbative series. Compar- 
ing to Eqs. (fT3)) and (fT4)) one also sees the sizable numer- 
ical cancellation between the two series at each order in 
as- The extent of this cancellation depends sensitively 
on the value of p'^"*. 

The perturbative uncertainties on these inclusive cross 
sections can now be used to determine the exclusive cross 
section uncertainties. Varying the scale up and down 
by a factor of 2 around mn /"i gives for the Tevatron 
CTtotai = (0.386 ± 0.040) pb and cr>i = (0.132 ± 0.034) pb 
with the pl^ and if'^'^ cuts as in Eq. (IT5|) . Adding these in 
quadrature according to the upper-left entry in Eq. (|4]) 
gives 



CTo = (0.254 ± 0.052) pb. 



(17) 



i.e., a 20% uncertainty. In contrast, when doing a 
scale variation directly in the fixed-order expansion for 
cro(p'^"*), as in Eq. ([1]), one implicitly assumes that the 
perturbative uncertainties between the series for crtotai 
and cr>i are 100% correlated, giving cto = (0.254±0.006). 
Here this leads to an underestimate for the remaining un- 
certainty. For the LHC, using the reference cuts, we get 
CTtotai = (8.70±0.75) pb and a>i = (3.08±0.59) pb, lead- 
ing to 



CTo = (5.63 ± 0.96) pb, 



(18) 



i.e., a 17% uncertainty. In contrast, the direct scale vari- 
ation for (To yields ctq — (5.63 ± 0.15), which is again an 
underestimate. 

The two procedures of evaluating uncertainties can be 
compared as a function of p™*, and in the upper left 
panel of Fig. [T] we do so for cro(p™*) for Higgs produc- 
tion. Results for cto are obtained at NNLO for the LHC at 
Ecra = 7TeV, using MCFM to calculate the p^."' depen- 
dence, FEHiP '30, ,31] for the total NNLO cross section, 
and /i — ran 12 for central values. The central value is the 
solid blue curve, and the green dashed and dotted lines 
show the results of direct exclusive scale variation by a 
factor of 2. For small values of p™' the cancellations that 
take place for (To(p'^"') cause the error bands to shrink. 
In particular, the direct exclusive scale uncertainty van- 
ishes at p^"' ~ 25 GeV, where there is an almost exact 
cancellation between the two series in Eq. (ITU)) , and the 
uncertainty curves pinch together. In contrast, the outer 
red solid lines show the result of our method, which com- 
bines the independent inclusive uncertainties to obtain 
the exclusive uncertainty, Ag = A^^^g^j + A?,j^. One can 
see that for large values of p^"' this combined inclusive 
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FIG. 1: Perturbative predictions for H + jets (upper left panel), WVK + O jets (lower left panel), H + 1 ]et with p'^^ > 30 GeV 
(upper right panel), and -ff + 1 jet with p^J > 120 GeV (lower right panel). Central values are shown by the blue solid curves, 
direct scale variation in the exclusive jet bin by the green dashed and dotted curves, and the result of combining independent 
inclusive uncertainties to get the jet-bin uncertainty by the outer red solid curves. 



uncertainty estimate reproduces the direct exclusive scale 
variation, since cr>i(p'^"') becomes small. On the other 
hand, for small values of p™' the uncertainties obtained 
in this way are now more realistic, because they explicitly 
take into account the large logarithmic corrections. The 
features of this plot are quite generic. In particular, the 
same pattern of uncertainties is observed for the Teva- 
tron, when we take /i = uih as our central curve with 
/i = 2mH and /z = to///2 for the range of scale variation, 
and whether or not we only look at jets at central rapidi- 
ties. We also note that using independent variations for 
Hf and /ir does not change this picture, in particular the 
Hf variation for fixed fir is quite small. 

Since both NLO and NNLO resuhs for cro(PT*) are 
available, it is also useful to consider the convergence, 
which we show in Fig. [2] for the Tevatron (top panels) 
and the LHC at 7 TeV (bottom panels) . In the left pan- 
els we directly vary the scales in (To{pf^^) to estimate the 
uncertainty, while in the right panels we again propagate 
the uncertainties from the inclusive cross sections. As we 
lower p™' , the direct exclusive scale variation uncertainty 



estimate decreases at both NLO and NNLO, and eventu- 
ally becomes very small when the curves pinch and the 
uncertainty is clearly underestimated. In contrast, the 
combined inclusive scale variation gives realistic uncer- 
tainties for all values of p™*. In particular, there is con- 
siderable uncertainty for small p™* where the summation 
of logarithms is important. 



B. Higgs + 1 Jet 



As our next example we consider the 1-jet bin in Higgs 
production from gluon fusion. This jet bin is defined by 
two cuts, one which ensures that the jet with the largest 
Pt is outside the 0-jet bin, pl^l > p™/, and one which 
ensures that the jet with the next largest px is restricted, 
Pt2 — Pt"'' ^o ttiat we do not have two or more jets. The 
1-jet cross section can be computed as a difference of 
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FIG. 2: Fixed-order perturbative uncertainties for gg — >■ _H"-|-0 jets at NLO and NNLO. The upper panels are for the Tevatron 
and the lower panels for the LHC with Ecm = 7TeV. On the left, the uncertainties are obtained from the direct scale variation 
in (To(p7^*) between /i — mH/4 and n = mn- On the right, the uncertainties are obtained by independently evaluating the 
scale uncertainties in crtotai and (T>i(p'^"*) and combining them in quadrature. (For the LHC case the dark shaded NNLO bands 
correspond to results in the top left panel of Fig. [l] The direct exclusive scale variation band corresponds to the dashed green 
lines, and the combined inclusive uncertainty band corresponds to the solid red lines.) 



inclusive cross sections with these cuts, 

ai = a>i(pjp^* > p^r\') - ^^Ml^Ti > PtUtI > Pt') • 

(19) 
For convenience we adopt the notation that p™* is always 
used for the cutoff that determines the upper boundary 
of the jet bin under consideration, which gives the analog 
of the L dependent terms in Eq. ((TU)) . 

The inclusive cross section a>i that includes the 1- 
jet bin exhibits large perturbative corrections, much as 
Ctotai does for the 0-jet bin. For cr>i the large corrections 
are caused in part by the large double logarithmic series 
in hilpi^^/mH), but remains predominantly independent 
of the large double logarithms oi L = ln(p')^2/"^-ff ) which 
control the series for cr>2. With jif = fir = ruH /'2, rnn = 
165 GeV, and MSTW2008 NNLO PDFs, we find 



Jet 



cr>i(p^Ti > 30 GeV) 



o->2(Pn > 30GeV,pJ^2 > 30 GeV) 
= (2.00 pb) [3.6 as + 0(a^)] . 



(20) 



(2.00 pb) [1 + 5.4 0,-1-0(02 



For (Ti = (T>i — (7>2 there is a sizable cancellation between 
these Us terms. If we lower the cut to p!^2 ^ 22 GeV then 
the logarithm increases and there is an almost exact can- 
cellation with the 5.4as. In the top right panel of Fig. [T] 
we plot (Ti as a function of p™*, and we again see that this 
cancellation occurs in a region where there is a dramatic 
decrease in the direct exclusive scale dependence (green 
dashed and dotted curves). Using the inclusive uncer- 
tainties for (7>i and (T>2, and adding them in quadrature, 
gives the outer solid red curves, which again avoids this 
problem and provides a more realistic estimate for the 
perturbative uncertainty. 

Using the result from the Appendix we can examine 
the full uncertainties and correlation matrix with 0, 1, 
and (> 2)-jet bins in Higgs production. For the cuts 
in Eq. (PU)) varying the scale by factors of 2, we have 



CTtotai = (8.70 ± 0.75) pb, cr>i = (3.29 ± 0.62) pb, and 
('■>2 — (0.85 ± 0.49) pb, corresponding to relative uncer- 
tainties of 8.6%, 18.8%, and 57%, respectively. We let 
b(x) denote the relative percent uncertainty of the quan- 
tity X, and p[x, y) the correlation coefficient between x 
and y. The Appendix yields 



(0.31 pb) [3.7a,, ^-©(a^)] , 



(24) 



,5(cto) = 18%, 

p(fT0,O'total) = 0.77, 

p(cro,o-i) = -0.50, 



<5(ai) = 32%, 
piai,a>2) = -0.62, 



(21) 



where we have only shown the nonzero correlations. Note 
that <7q and txi as well as cti and cr>2 have a substantial 
negative correlation because of the jet-bin boundary they 
share, while (Tq and a>2 are uncorrelated. 

In contrast, the direct exclusive scale variation results 
in all the cross sections being 100% correlated. Be- 
cause of the cancellations between the perturbative se- 
ries, this leads to much smaller (and unrealistic) un- 
certainties, with our choice of cuts S{(Jo) — 2.3% and 
S{ai) = 5.5%, which is reflected in the pinching of the 
dotted and dashed green lines in Fig. [T] (Note that in- 
creasing the range of scale variation or separately varying 
fj,r and II f does not mitigate this problem.) The analog 
of Eq. ^ for this example would be 

0.62 X 2.3% -I- 0.28 x 5.5% -I- 0.10 x 57% = 8.6% . (22) 

When all ai are 100% correlated, ctq is forced to have 
a smaller relative uncertainty than crtotai, as in Eq. (jlj, 
since it has to make up for the much larger uncertainties 
in (j>2. 

In addition to the cross sections in each jet bin, we can 
also consider the relative jet fractions /o — cro/atotai and 
fi/ctotai, which are often used in experimental analyses. 
The perturbative theory uncertainties and correlations 
for the jet fractions follow by standard error propagation 
from those in Eq. (I21[) . The general expressions are given 
in the Appendix, and we find 



^(/o) = 13%, 
p(/o,o-totai) =0.42, 
p(/o,/i) = -0.80. 



<5(/i) = 33%, 

P(/l,crtotal) = -0.26, 



(23) 



Comparing to Eq. (|2T|1 . the use of jet fractions with fXtotai 
in the denominator yields a nonzero anticorrelation for 
ctotai with the 1-jet bin, and decreases the correlation 
for CTtotai with the 0-jet bin. 

It is also interesting to consider the case with p!^^ > 
120 GeV, where the logarithms oi pl^j^/mn are not large. 
The cross section a>i now has a smaller perturbative 
correction, but for a region of cuts on p!^2 there are still 
substantial cancellations in cti. For instance, for p!^^ ^ 
60 GeV we have 

cr>i(p^Ti > 120 GeV) 

= (0.31 pb) [1 + 2.9 a, + 0{al)] , 

cr>2(Pn > 120GeV,pJ^2 > 60 GeV) 



and the a^ terms completely cancel around p!^2 — 
70 GeV. In the bottom right panel of Fig. [T] we plot 
(7i as a function of p™* for this scenario. Once again the 
combined inclusive uncertainties (solid red curves) give 
a better estimate than the direct exclusive scale uncer- 
tainty determined by up/down /i variation in cti (green 
dotted and dashed curves). It is interesting to notice that 
the curves dive and a logarithmic summation in p!^^ be- 
comes important earlier now, i.e., at much larger values 
ioi p!^2^ when the cut onp^^-^ is raised. For p!^^ > 120 GcV 
and pi^2 ^ 30 GeV fixed-order perturbation theory does 
not yield a controlled expansion, and the resummation of 
the jet- veto logarithms is clearly necessary. 



C. WW + Jets 

The process pp —> WW + jets is the dominant irre- 
ducible background for the H — )■ WW* search in the 0-jet 
bin, and also exhibits a relatively large K factor ~ 1.5. 
Hence, it is interesting to contrast the scale uncertain- 
ties here with those found ior H + jets. Including the 
Higgs search cuts (modulo the jet veto) , the K factor for 
WW becomes larger than 2 [l^l , but we will not include 
those cuts in our analysis here. With fir ^ fJ-f = fnw, 
NLO MSTW2008 PDFs, and a, = as{mw) = 0.1226, 
the total pp — )► WW cross section is 

atotai = (32.5 pb) [1 + 3.6 a, + 0{al)] , (25) 

while for the inclusive 1-jet cross section with logarithms 
of p™* we have 

<y>i{l^T > 30 GeV) = (32.5pb)[2.8a, + 0(a^)] . (26) 

Thus, when we consider ctq = crtotai — cr>i there is a siz- 
able cancellation for the a, terms. In Fig. [I] lower left 
panel, we show ctq for pp -^ WW + jets as a function 
of p™* . Once again the dotted and dashed green curves 
from direct exclusive scale variation exhibit a pinching 
near p™* ^ 30 GeV due to cancellations between the two 
perturbative series in Eqs. ([25)1 and ([25)1 . leading to an 
underestimate of the perturbative uncertainty. The com- 
bined inclusive uncertainty estimate again mitigates this 
problem. The pattern of uncertainties here is the same 
as for H + jets and -ff -I- 1 jet, just with smaller overall 
uncertainties. Just like for H + jets using independent 
variations for fif and fj,r does not change the picture, the 
fif variation for fixed fir is again quite small. 



D. W + Jets 

The exclusive process pp -^ W + N jets is an impor- 
tant benchmark process at the LHC and also an impor- 
tant standard model background for new physics searches 
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FIG. 3: Fixed-order perturbative uncertainties for the exclusive pp ^ W + 0,1,2 jet cross sections at NLO for the LHC with 
-Bern = 7TeV. Central values are shown by blue solid curves, direct exclusive scale variation in the exclusive jet bin by the 
green dashed and dotted curves, and the result of combining independent inclusive uncertainties to get the jet-bin uncertainty 
by the outer red solid curves. 



looking for missing energy. In this section we consider 
pp ^' W + jets, which provides us with a case to test 
our method when the perturbative corrections in the in- 
clusive cross sections are not as large. For simplicity, we 
only work to NLO here. Using fj,f ^ fir = "m-w for the 
central value and MSTW2008 NLO PDFs, the inclusive 
W production cross section is 

^totai = (80.7nb) [1 + 1.3a, + 0{al)] , (27) 

where we have summed over W , and have not included 
the leptonic branching fractions. For the inclusive 1-jet 
cross section we have 

o->i (^T* ^ 30 GeV) = (80.7 nb) [0.9 a, -1- ©(a^)] . (28) 



The perturbative coefficients in Eqs. (P7)) and (|28p are 
much smaller than in Higgs production. The resulting 
predictions for cro(Pr'') ^^^ shown in the top left panel 
of Fig. 131 where the different lines have the same mean- 
ing as in Fig. [TJ Since the as corrections are not very 
large here, the /i/ scale variation in the PDFs dominates 



over the /x^ variation in as and produces a 100% neg- 
ative correlation between crtotai and cr>i. (Keeping fif 
fixed at mw and only varying fir results in the expected 
pinching of the dotted and dashed green lines.) This 
means their scale uncertainties add linearly in ctq, which 
maximizes the uncertainty in this 0-jet cross section. In 
this case, our method, shown by the outer solid red lines, 
gives an uncertainty band very similar to direct exclusive 
scale variation. Hence, our method of using independent 
inclusive uncertainties still remains consistent for this sit- 
uation. 



E. W + 1 Jet 

For pp ^ W +1 jet the perturbative corrections in (T>i 
are larger than those in the W total cross section, which 
is in part influenced by logarithms from the lower cut on 
p!^^ , the pt of the leading jet. The situation for the W^-l- 1 
jet bin is similar to 7J -|- 1 jet. Considering Eq. ([T^ the 
series for the inclusive 2-jet cross section, cr>2, has large 



double logarithms L = ln(p^^2/'^w) of the second largest 
jet pt, which are independent of those in the perturbative 
series for cr>i. Taking /i — mw for central values, and 
using MSTW2008 PDFs at NLO, the total W+ + W 
cross sections with both jet cuts at 30 GeV are 

<7>i{l^Ti > 30 GeV) 

= (8.61 nb) [l + 3.4 a, + 0{al)] , 

<^>2{p't\ > 30GeV,pi^2 > 30 GeV) 
= (8.61 nb) [2.5 a^ + 0(02)] . 



(29) 



Once again the result for a>2 and the precise cancella- 
tion that occurs in cti is quite sensitive to p^^, the cut on 
p^^2; yielding an almost exact cancellation of the 3.4 a^ 
for pl^2 ^ 25 GeV. In the top-right panel of Fig. [3] we 
plot (Ti as a function of p^^ , with direct exclusive scale 
variation (green dashed and dotted curves) and those de- 
rived from independent inclusive uncertainties (solid red 
curves). Just like for H +1 jet, the direct exclusive scale 
variation curves pinch, while the inclusive curves avoid 
this problem and remain realistic. 

We can also consider what happens when we make a 
larger cut on p!^-^ . Here, unlike for the Higgs case, the rel- 
ative size of the perturbative correction in a>i increases. 
For instance, 

c^>i(p^Ti > 80 GeV) 

= (1.07 nb) [1 + 5.3 a, + 0(a^)] , 

c^>2(p^Ti > 80GeV,pJ^2 > 60 GeV) 

= {1.07 iih)[4:.l as + 0{al)]. (30) 

For p'^2 ^ Pt"* in- cr>2 the resulting 1-jet cross section cri 
is shown as a function of p^"' in the bottom-left panel of 
Fig. [31 The situation for the uncertainties is similar to 
that for the less stringent cut on p-!^-^ in the upper-right 
panel. Much like in 7J -I- 1 jet the logarithms start to 
influence the cross section at larger values of p^"* for the 
larger plp^ cut. 



F. W + 2 Jets 

As our last example we consider W + 2 jets, and for 
simplicity we only consider the case of W'^ production. 
The inclusive 2-jet and 3-jet cross sections with all jets 
cut at 30 GeV are 

cT>2{l^Ti.2 > 30 GeV) = (1.60 nb)[l + 1.0 a^ + C'(a^)] , 

^>3 (pti 2 3 > 30 GeV) = (1.60 nb) [2.3 a^ + 0(a^)] , 

(31) 

and the resulting exclusive 2-jet cross section as a func- 
tion of the p^* on the third jet is shown in the bottom- 
right panel in Fig. [31 

There are two different types of diagrams contribut- 
ing to this process, those having two external quark lines 



and two gluon lines at lowest order (qqgg), and those 
having four external quark lines at lowest order (qqqq). 
The qqgg-type contributions have the same behavior as 
PF -|- 1 jet, again displaying a pinching in the direct ex- 
clusive scale variation curves. On the other hand, in 
the qqqq-type contributions the PDF scale dependence 
dominates, similar to what we observed for VF -|- jets. 
The combination of the two leads to the behavior seen in 
Fig. [3lat large Pj^^, where the scale uncertainties in the 
inclusive 2-jet cross section are asymmetric. Here there 
is some choice for how to combine the scale variation into 
an uncertainty estimate for a>2 (green dashed and dot- 
ted curves). The choice one makes for a>2 simply propa- 
gates into the equivalent choice for the exclusive 2-jet bin 
a2 (solid red curves). For simplicity in the bottom-right 
panel of Fig. [31 we still use fi = mw/2 and /i = 2mw to 
determine A>2, in which case the central value should 
be taken as the center of the band rather than the dark 
solid blue line for /i = mw- 

For W+2 jets in Fig.[31the pinching caused by the qqgg 
contributions is again mitigated by combining the inclu- 
sive uncertainties. Hence, we see that our method can 
be applied and gives more stable uncertainty estimates 
even in more complicated cases where several components 
contribute to the cross section. 

Note that we have also checked that when increasing 
the cuts on the two leading jets, the same effect as in H+1 
jets and W^ 4- 1 jets happens here as well. Namely, the 
jet- veto logarithms from restricting the third jet become 
more important earlier and influence the cross section at 
larger values of p™* for larger pi^^ ^ cuts. 



IV. RESUMMATION FOR HIGGS + JETS 

In Sec. mil we have seen that direct exclusive scale 
variation often leads to an accidental underestimate of 
the uncertainties for exclusive jet bin cross sections for 
a range of experimentally relevant cuts. Instead combin- 
ing independent uncertainties on inclusive cross sections 
yields a more uniform (and larger) uncertainty band for 
the exclusive jet bins. The region where direct exclusive 
scale variation runs into trouble borders the region where 
the resummation of the large logarithms of •p'!^ becomes 
important. In this section, we test how realistic the fixed- 
order scale uncertainties are by comparing them to a case 
where the resummation of large logarithms induced by 
the jet bin are known to NNLL-(-NNLO accuracy. 

We again consider 7f -I- jets from gluon fusion. At 
NNLL order accuracy the resummation is sensitive to 
the precise jet algorithm used to define p'!^ , and other 
complications in the required theoretical setup. To avoid 
these issues, we will use a slightly different variable to 
define the 0-jet bin, an inclusive event shape known as 
beam thrust kMII. 



rcn. = ^(i?fc-b^| 



(32) 
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FIG. 4: Comparison oi gg ^ H + jets using T^^"' at fixed NNLO with the resummed results at NNLL+NNLO. For the 
fixed-order uncertainties in ctq we use independent inclusive scale variations in atotai and ct>i. The uncertainty method for the 
resummed results is described in the text. 



The sum over k runs over all particles except the Higgs 
decay products. Beam thrust essentially measures the 
thrust of an event along the z beam axis. When T^^j < 
7^^*, from Eq. (|32l) we see that events in cro(7^m) are only 
allowed to contain hard radiation in the forward regions 
at large rapidities, and hence this cut vetoes central jets. 
Much like with p^ the perturbative series for this aa has 
double logarithms, for example the analog of Eq. ^ is 



•in, T-cut 
„ /T-cutN „ I 1 "^"s 1 2 'cm 
C^OI /cm ) = CTB 1 In 



+ ■ 



(33) 



For beam thrust, the all-order resummation of per- 
turbative corrections is known to NNLL order for both 
H + Q jets and ^ -I- jets [H-H^]. For Higgs produc- 
tion the computation has been extended to fully include 
all NNLO corrections, and it was observed that the re- 
summed cross section at NNLL-I-NNLO had larger uncer- 
tainties than the pure NNLO result for crolTcm < Tcm) 
utilizing direct exclusive fixed-order scale variation. This 
led to the conclusion that the direct exclusive scale vari- 
ation underestimates the fixed-order perturbative uncer- 
tainties in the 0-jet bin. In the resummed calculation, 
fixed-order as expansions are carried out at three dis- 
tinct scales (hard ^h, jet /beam ^b, and soft ^s) which 
appear in the corresponding factorization theorem. The 
uncertainties in the resummed cross section are obtained 
by varying these scales. Varying /i/f up and down by a 
factor of 2 moves all three scales up and down, and hence 
is a scale variation that is correlated with the usual scale 
variation for the inclusive cross section. Varying /i^ or 
/i5 while holding ^h fixed explicitly accounts for addi- 
tional higher order uncertainties induced by the presence 
of the large jet-veto logarithms, and hence allow us to 
determine Acut- 

In Fig. m we compare the remaining perturbative un- 
certainties after resummation at NNLL-I-NNLO, shown 
by the darker red bands, to the NNLO uncertainties 



obtained with the fixed-order method advocated here, 
which are shown by the lighter gray bands. The results 
for the NNLL+NNLO cross section are obtained from 
Ref. [33|.^ The left panel shows the results for the Teva- 
tron and the right panel the results for the LHC at 7 TeV. 
The fact that the resummation reduces the perturbative 
uncertainties, as it should, shows that our method of us- 
ing independent inclusive scale variations yields more ro- 
bust fixed-order uncertainties. 



In the resummed calculation, Ctotai is by construction 
not affected by the /ig and ^s variations. We denote 
the combined ^s and ^b uncertainty by A5B. It pro- 
vides a direct estimate of the cut-induced uncertainty, 
Acut = A5B, which is anticorrelated between croCT^m*) 



■cut^ 
cm 



and the corresponding cr>i(7^™*) = crtotai - cro(7^' 
On the other hand, the ^h variation affects all the cross 
sections yielding an uncertainty component that is 100% 
correlated between them. In particular, it is responsi- 
ble for estimating the perturbative uncertainty of CTtotai, 
for which it is equivalent to the usual fixed-order scale 
variation, l^utot = Atotai- The full covariance matrix for 
{ crtotai, ctq J o->i}, that is the analog of Eq. p^ but for 



We have made a small improvement to Ref. 13711 . The 
NNLL+NNLO results of Ref. [H fuUy incorporate the NNLO 
corrections by adding so-called nonsingular fixed-order contribu- 
tions, which are terms that do not appear in an expansion of the 
strict NNLL result. In Ref. |37| the nonsingular contributions 
were obtained for the sum of 0(as) + 0{a^) cross-sections using 
FEHiP |30| . |31| . Here we use a much higher statistics spectrum 
from MCFM [17| . which allows us to separately determine the 
nonsingular cross sections at 0{p.s) and 0{al). The only place 
this improvement is visible is for T^'jJJ* < 3 GeV, where the re- 
summed cross sections are now consistent with zero within the 
displayed uncertainties. 
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the resummed result, is then 
C = CsB + Ch , 



CsB ~ 






\ 




Al^ 


'-^SB ' 




-A|5 


A|b/ 




^Ito, 


A//totAHo 


A//tot Ah>i 


A//tot A//0 


AL 


Ai/o Ah>i 


A/ftot Ah>i 


A_ffo Ah>i 


Al,>i 



(34) 



Cff = 



where IS.sb is obtained from the envelope of the //g and 
/is variations, and 6*55 is equivalent to Ccut in Eq. ([5]). 
The A/fi are obtained from the /z/f variation and satisfy 
^Htot = A/fo + A/f>i. The full uncertainty in the 0-jet 
bin shown by the darker red bands in Fig.|4]is then given 
by A|g + A|j^g , which is the 0-bin entry on the diagonal 
ofC.6 

Compared to Eq. (|M|) . using a direct exclusive scale 
variation at fixed order would correspond to taking 
A5B -^ and obtaining the analog of the lS.Hi by scale 
variation without resummation {puj = fJ-B = Ms)- On the 
other hand, our proposed fixed-order method would cor- 
respond to taking Asb -^ A>i and Ah>i -^ 0, such that 
Ahq — A/ftot -^ Atotai- Hence, the resummation of the 
jet-veto logarithms allows one to capture both types of 
uncertainties appearing in the two different fixed-order 
methods. Note that the numerical dominance of A|^ 
over AhqAh>i in the 0-jet region is another way to jus- 
tify the preference for using the combined inclusive scale 
variation over the direct exclusive scale variation when 
given a choice between these two methods. 

As an example, consider 7^™' = 20GeV. At fixed 
NNLO, the inclusive cross sections are crtotai — (8.70 ± 
0.75) pb and (t>i = (2.25 ± 0.62) pb. Using Eq. ^, this 
gives 



<5(cro) = 15% , 
/o(o-o,crtotai) = 0.77, 
p{ao,a>i) = -0.64. 



J(a>i) = 28%, 

p(a>i,CTtotal) = , 



(35) 



For (To this corresponds to the lighter gray bands in Fig. 21 
and the structure here is very similar to what we saw in 
Eq. dUD- 

From our resummed result using Eq. p4p we obtain 

Siao) = 11.8% , S{cT>i) = 19.7% , 

P(a'0, CTtotal) = 0.04 , p(o->l, CTtotal) = 0.33 , 

/o(^o,CT>i)--0.82, (36) 



In the results of Ref. (37j , the envelope of all three scale variations 
was used to obtain the total uncertainty. The slightly modified 
procedure we use here, which adds Asb and Ah in quadrature, 
gives very similar results, but has the advantage that it also 
allows for a straightforward treatment of the correlations. 



which for ctq corresponds to the darker red bands in 
Fig. 21 After resummation neither of uo and a>i is 
strongly correlated with (Ttotai anymore, which at first 
sight is perhaps a bit surprising. However, for small 7^*^"' 
this is not unexpected and is simply due to the fact that 
the central values and remaining perturbative uncertain- 
ties are dominated by the resummed logarithmic series 
(i.e. Asb dominates numerically over Aho and Ah>i)- 
In fact, this supports our arguments in Sec. [TTl that the 
uncertainties from higher-order terms in the logarithmic 
series for cr>i and the fixed-order series for atotai can and 
should be considered independent, which lead to Eq. (|lip . 
Comparing Eqs. ((55)) and ([55)) . we see that the uncer- 
tainties obtained from our fixed-order method follow a 
similar pattern for the relative uncertainties for ctq and 
(T>i as observed in the resummed result, with a strong 
negative correlation between them. Since resummation 
provides an improved treatment of the cut-induced ef- 
fects, we take this as further evidence that the method 
of using inclusive fixed-order cross section uncertainties 
provides a consistent way to obtain reliable estimates of 
perturbative uncertainties in exclusive jet bins. In par- 
ticular it provides a suitable starting point for an un- 
certainty estimate, that can be further refined when an 
appropriate resummed result becomes available. 



V. CONCLUSIONS 

We have proposed a method to estimate perturbative 
uncertainties in fixed-order predictions of exclusive jet 
cross sections that accounts for the presence of large 
logarithms at higher orders caused by the jet binning. 
The method uses the fixed-order calculations of inclusive 
cross sections, (t>n and a>]\[^i, for which the standard 
scale variation provides reasonable uncertainty estimates, 
and combines these inclusive uncertainties into an esti- 
mate for the corresponding exclusive A'^-jet cross section 
(Tat = a>N — o'>N+i, treating the inclusive cross sections 
as uncorrelated. 

We have illustrated this procedure for a variety of pro- 
cesses, including analysis of iJ -t- 0, 1 jets, WW + jets, 
and W + 0,1,2 jets with MCFM, and showed that it 
yields more robust estimates of theory uncertainties than 
direct exclusive scale variation. We have also shown for 
a specific case with H + jets that it leads to fixed-order 
uncertainties that are theoretically consistent with the 
corresponding resummed predictions. In jet bins used 
for new physics searches, we anticipate that it should 
yield realistic uncertainty estimates for standard model 
backgrounds. We also expect that it provides a suitable 
fixed-order starting point for the central values, uncer- 
tainties, and jet bin correlations, which can be improved 
by higher-order logarithmic resummation. 

Our treatment of the fixed-order exclusive and inclu- 
sive cross sections has followed the standard approach of 
always using cross section results at the same order in 
as- It would be interesting to study whether this can be 
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relaxed when using differences of inclusive cross sections 
to compute the central values for the jet bins. For exam- 
ple, for gg ^ H one could independently compute fXtotai 
at NNLO, and cr>i and a>2 each at NLO, and then use 
these to compute the jet bins as cto ~ crtotai — cr>i and 
CTi = cr>i — cr>2. Since we argued that the inclusive se- 
ries can be treated independently, it may be consistent to 
include them to different orders to compute the central 
value and uncertainties of ui . This would have the ad- 
vantage of allowing one to utilize the NLO result for cr>2 
without destroying the consistent perturbative expansion 
for CT>i and atotai when the jet bins are added together. 
Since in this case the perturbative order of the jet bound- 
ary between cti and ct>2 does not match up, this deserves 
a dedicated study before being used in practice. 
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Appendix: Case of Three Jet Bins 



In this appendix we generalize Eq. p2|) to the case of 
0, 1, and (> 2) -jet bins that is actually used in current 
Higgs searches. Since only neighboring jet bins will be 
correlated, the generalization to more than three jet bins 
is not any more complicated. 

We start from the inclusive cross sections iJtotai, cr>i, 
cr>2, and denote their absolute uncertainties by Atotai, 
A>i, A>2 and their relative uncertainties by 5i = lS.i/ai. 
We define the exclusive cross sections and event fractions 
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The covariance matrix for 
{o-totai, cfq, 0-1, cr>2} is given by 
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The relative uncertainties for /o and /i are 

5{hf^{j^-l)\5l,a, + 5l,), (A.5) 

5{hf = Cai + (^)'^>i + (^ - l)'42 , 

and their correlations are 
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